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The domain wall fermion chiral condensate in quenched QCD 
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We examine the chiral limit of domain wall fermions in quenched QCD. One expects that in a quenched 
simulation, exact fermion zero modes will give a divergent, 1/m behavior in the chiral condensate for sufficiently 
small valence quark masses. Unlike other fermion formulations, domain wall fermions clearly demonstrate this 
behavior. 



1. INTRODUCTION 

It is a common assertion that quenched lat- 
tice QCD qualitatively produces many of the fea- 
tures found in continuum full QCD. However, a 
quenched simulation is expected to include topo- 
logical gauge configurations which produce exact 
fermion zero modes. Not suppressed by the miss- 
ing fermion determinant, such modes will cause 
divergent behavior in the small mass limit of the 
chiral condensate. 

Surprisingly, the standard lattice fermion for- 
mulations do not demonstrate divergent behavior 
for the quenched chiral condensate at current lat- 
tice spacings and quark masses. One may imagine 
that the small eigenvalue spectrum of those lat- 
tice Dirac operators bears little resemblance to 
the continuum Dirac spectrum until the lattice 
spacing is substantially smaller. See the introduc- 
tion of H] for further discussion and references. 

Starting from Kaplan's initial ideas for domain 
wall fermions (DWF) [|, several fermion formu- 
lations have been developed which use an infi- 
nite number of massive fermion flavors to pro- 
duce an effective chiral-invariant action for a sin- 
gle massless fermion flavor. Using the overlap 
formalism Q], which describes the infinite flavor 
case, Narayanan and Neuberger demonstrated 
that these theories have a unique, integer valued 
index for the associated Dirac operator. Stud- 
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ies have shown that nearly exact zero modes per- 
sist in truncated theories, provided the number 
of heavy flavors is large enough. One should then 
expect to see these zero mode effects in quenched 
simulations. 

2. THE CONTINUUM PICTURE 

Using a spectral decomposition of the contin- 
uum Dirac operator, we can express the quark 
chiral condensate as a function of the quark mass 
for any ensemble of gauge fields 
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where we have used the 75 symmetry of the Dirac 
spectrum and N zm are the number of zero eigen- 
values in the spectrum. We can similarly express 
the integrated pion correlator 
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Spontaneous chiral symmetry breaking is related 
to a non-zero eigenvalue density at A = by the 
Banks-Casher relation which, for the case of non- 
zero N zm , becomes 
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+ ?rp(0) +0(m), m->0. (3) 



Finally, we note that the identity 
m [ d i x(Tr(x)ir{0)) = (^) 



(4) 



follows from the continuum chiral symmetry as 
represented by equations (P and 

3. DOMAIN WALL FERMIONS 

Our numerical simulations were performed us- 
ing the boundary fermion variant of domain wall 
fermions ||. Here one adds a fifth dimension to 
the normal 4-D space time of extent L s . The 
gauge fields are not changed by this extension 
and do not depend on s. There are also two 
mass parameters for the fermions: mr, the ex- 
plicit four-dimensional bare quark mass and m , 
the five-dimensional bare quark mass. Except for 
different symbols for the three parameters listed 
above, our conventions follow ||. 

In the free field case, the effective mass of the 
low energy state is ||] : 

mfg cc) = m (2 - m ) [m f + (1 - m ) is ] . (5) 

where, for one flavor physics, mo should be in the 
range (0, 2). This simple equation reveals two im- 
portant features of domain wall fermions. First, 
as L s — > oo, m c ff a m/ and the proportionality 
constant depends on the five-dimensional mass 
mo. Second, for finite L s , the Wilson-like addi- 
tive correction to m e g is exponentially suppressed 
in L s and the rate of decay is determined again 
by itiq. In the interacting case, we expect that 
the effective mass will behave in a similar fashion 
where the residual mass term m ros should vanish 
exponentially in L s at a rate governed by mo and 
the coupling constant. Indeed, numerical studies 
Q indicate that the decay remains exponential 
with a decay rate that becomes faster as the con- 
tinuum limit is approached. 

4. QUENCHED DWF AT ZERO TEM- 
PERATURE 

We begin by examining the extent to which do- 
main wall fermions reproduce the effects of chiral 
symmetry for finite L s in a quenched simulation 
by computing the ratio of the two sides of (Q) on 
8 3 x 32 lattices at (3 = 5.7 and 5.85. We choose to 
work at a domain wall height, mo = 1.65. This 
choice is supported by numerical studies of the 
chiral condensate for several values of mo, L s and 
mf at these values of /3. 
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Figure 1. Ratio R of m J d 4 x{Tr(x)n(0)) to (jpip) 
In Figure 1 we show the ratio of quantities pro- 
portional to the left and right hand sides of (^). 
The constant, m/-independent, behavior of this 
ratio, required by chiral symmetry, is increasingly 
visible in Figure 1, for the curves with larger L s . 
Particularly striking is the comparison of this ra- 
tio between (3 = 5.7 and 5.85. The weaker cou- 
pling with L s = 16 shows a degree of chiral sym- 
metry present at the stronger coupling only for a 
three times larger L s . This is consistent with the 
picture that stronger coupling requires larger L s 



quenched 8 3 x32, p=5.85, 
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Figure 2. 1/m/ divergence in (tptp) 

Next we study the small mass limit of (i/ji[>) 
computed on 200 quenched configurations at f3 = 
5.85 and L s = 32. The results, shown in Figure 
2, show a clear signal for a 1/m/ divergence for 
mf < 0.01. The solid line is a fit for 0.001 < 
mf < 0.1 which includes a 1/m/ term and the 
dashed line is a linear fit for 0.01 < m/ < 0.1. 
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The constant and linear terms are consistent be- 
tween these two fits within errors and the x 2 /dof 
is one for both fits. Table 1 summarizes the rest 
of our T — results. We have used mo = 1.65 for 
all simulations. 

We believe that this is the first time that the 
1/m/ behavior in expected from fermion 

zero modes present in quenched calculations has 
been seen. Notice the decrease in the strength of 
this term, shown in Table 1, when the volume is 
increased from 8 3 x 32 to 16 3 x 32. 



Table 1 _ 
Fit to (ipip) 



-i/mf + Co + ci? 



8 3 x 32 
(3 = 5.70 
L s = 48 
40 configs 



8 3 x 32 
(3 = 5.85 
L s = 32 
200 configs 



16 3 x 32 
(3 = 5.85 
L s = 32 
90 configs 
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quenched 16 x4, m =1.90, L s =32 



c_i(xl0" e ) 


1.6(4) 


3.8(3) 


0.60(9) 


co(xl0- 3 ) 


2.15(3) 


0.82(3) 


0.89(1) 


ci(xl0- 2 ) 
X 2 /dof 


6.01(2) 


6.73(3) 


6.71(2) 


8.99 


1.24 


2.96 



5. QUENCHED DWF AT FINITE TEM- 
PERATURE 

In Table 2 we summarize a quenched study of 
(4>ip) on a 16 3 x 4 lattice evaluated for (3 slightly 
above (3 C = 5.6925. Note, we continue to see a di- 
vergent behavior of (iptp) ■ To our surprize, we also 
see a non-zero constant contribution to the chiral 
condensate well into the symmetric phase. Figure 
3 shows the behavior of both c_i and Co as we in- 
crease (3 above (3 C . Each shows the rapid decrease 
as the system heats up expected if their origin 
is topology-induced fermionic, zero modes. This 
never-before-seen, quenched small-mass behavior 
of (ijjip}, suggests that domain wall fermions of- 
fer new insights into the continuum behavior of 
quenched QCD. 
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Figure 3. Fit coefficients for T > T c 
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